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Abstract 

This paper presents a methodology of detecting rotor imbalances, such as mass imbalance and crack-induced 
imbalance, using shaft synchronous vibrations. An iterative scheme is developed to identify parameters from 
measured synchronous vibration data. A detection system is integrated by using state-of-the-art commercial analysis 
equipment. A laboratory rotor test rig is used to verify the system integration and algorithm validation. A real engine 
test has been carried out and the results are reported. 
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complex symbol, j = V^I 
crack induced imbalance constant 
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Laplace transform of y(t). L( /oj) = L[ v(7)] 

time domain displacement response 
crack induced imbalance, U c =K c or 
mass imbalance in frequency domain 
amplitude of the imbalance, \U m \ 
phase angle of the imbalance 


damping ratio of a single-degree-of-freedom system the system, 
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shaft speed 

natural frequency of a single-degree-of-freedom system, 4k / m 


I. Introduction 

In rotating machinery, it has been well established that mass-induced (ref. 1) synchronous vibration is 
proportional to the square of the rotational speed, i.e., co 2 , in the region where a>«a>o, where a> 0 is the Gi st resonant 
frequency of the machinery system. It has also been established that radial-axial cracks (refs. 2 and 3) in the rotating 
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disk or shaft can cause a synchronous vibration that is proportional to co 4 , under the same assumption of o>«(Ao- 
Several groups have developed detection systems based on these assumptions (refs. 4 and 5). Unfortunately, in 
reality, there are several major difficulties in implementing these traditional algorithms: 


1 . It is difficult to satisfy the k>«c»o condition without significantly reducing the signal-to-noise-ratio, because 
the first natural frequency in many rotating machinery systems is very low. In soft-mounted rotor systems, the 
first natural frequency can be so low that it prevents such monitoring systems from getting meaningful 
results. 

2. It is difficult to simplify real rotating machinery into a single degree of freedom (SDOF) system. 

Additionally, the vibration sensors are usually mounted on the bearing case or the engine case, so that 
additional structural resonances, besides rotor-shaft resonances, may be picked up, thus complicating the 
simple polynomial relationship. 

In this research, a new algorithm is developed for mass imbalance and crack detection in a multi-degree of 
freedom (MDOF) rotating machinery system. The proposed method takes advantage of high signal-to-noise-ratio in 
the neighborhood of resonant frequencies in synchronous vibration, thus it is expected to have better accuracy in 
imbalance detection. 


In the simplest possible model for rotor vibration analysis, the system is described by a spring-mass-damper 
system as shown in figure 1. With this model, the system has to be either a rigid rotor and disk with a flexible 
bearing, or a flexible mass-less rotor shaft with a rigid bearing and single rigid disk. Under these assumptions, the 
frequency domain displacement response due to a mass imbalance can be expressed as 


The presence of certain cracks in the rotor disk introduce similar vibration to that caused by a mass imbalance, 
i.e., 


II. Theory 



(i) 



C/ c ro 2 


( 2 ) 



C B = bearing damping 

C s = shaft damping 

d e = distance 

E = shaft modulus 

I = shaft moment of inertia 

K b = bearing stiffness 

M e = unbalance mass 

m = rotor mass 

u m = unbalance 

9 = angular position 

to = speed 



Rigid Rotor, Flexible Bearing 


Mass Unbalance 


Rigid Bearing, Flexible Rotor 


Single-degree-freedom Dynamics Model 


Figure 1. — SDOF rotor dynamics modeling. 
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where Uc is the crack induced imbalance. It is a complex variable in the frequency domain. According to reference 3 
and finite element verification in reference 6, Uc is a function of rotational speed, crack length, and rotor disk 
geometry, or more precisely. 


U c =K c ar =ifc(crack, disk)® 2 (3) 

Thus, in presence of both mass and crack- induced imbalances, the rotating system possesses a displacement 
vibration response format as 


Y(M=- 


U m (o 2 


K C (0 4 


/n (©0 - © 2 + 2 yl;o©o©) w(©q - © 2 + 2 y'l;o©o©) 

In the region where ©«© 0 , equation (4) can be simplified as 
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are the imbalance and system parameter related complex constants. For a rotating system operated in a speed range 
between resonant frequencies, an improved detection methodology can be found in reference 7. 

Real rotating machinery can rarely be modeled as a simple system. It is usually a distributed system with 
infinitely many DOFs and with low first resonance frequencies. In such cases, the system response has to be 
modified into a multi-degree-of- freedom (MDOF) format. Assuming the system has distinct natural frequencies, for 
simplicity, we have the displacement response as 


Y(M=Y, 

i = 1 
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Provided that the rotor is operated in the frequency range a> = [oo 0 i, ram], and 

K>oi < (0< ©02 

©a- < ®oi < ©a+1 and ©a +/ < ©02 < ©A+/+1 

k> 1 or ©a = 0 

/>! 


and further supposing that all the resonant frequencies are well separated such that the cross coupling from modes 
outside the operational range are in a simple form, then equation (7) can be expressed into a finite summation as 


>"(./©) = x 


f7.ii/© 


/ 


Kaco 4 


i=k 


»;, (©? -© 2 + 2 y '^©/©) ^ mj (©? - © 2 + 2 y '^©/©) 


+ C 


(9) 


where C is a complex constant to account for the influence from other modes. For a given response curve, ©,- and c„ 
i=k, k+1, .... I, can be approximated by identifying the peak location and by the half power method to get initial 
estimations. Then an iterative method, such as the diagram shown in figure 2, can be implemented to improve the 

JJ 

accuracy and to solve for relative modal excitations of — — , and — — 

/«, m i 
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Figure 2. — Diagram of solving imbalances. 


In cases when accelerometers are used as sensors, equation (9) needs to be modified as 
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(10) 


where C is a complex constant representing the influence from other modes, and Z),s are used to account for 
responses due to random excitations imposed to the rotating system, such as environmental excitations and nonlinear 
interactions of the system. 


III. System Integration 

In an accelerated test, the engine ramp-up (or coast-down) process is actually taking place in a very short time. 
It is very important to have a system that can handle multi-channel simultaneous data acquisition, synchronous data 
processing, condition monitoring, and data management in a short time. Additionally, the signals from such tests are 
usually noisy, thus it is critical to have a processing technique which can improve the signal-to-noise-ratio. 

To meet the engine requirements, the Pulse system by B&K was selected as the hardware platform. The 
Production Test Advisor (PTA) by Signalysis, Inc. was selected as the database management system. A MATLAB 
(The MathWorks, Inc.) based post-processing system was developed in-house at GE and integrated into the system 
for real-time health monitoring. A brief system diagram is shown in figure 3. 

Sensors 

Sensors include vibration transducers and a tachometer signal. For synchronous vibration detection, any 
vibration sensors, such as accelerometers, velocity sensors, or displacement probes, can be used. Appropriate signal 
conditioners are needed for the corresponding sensors used. A once per revolution (1/rev) signal from the rotating 
shaft is also needed for analysis and order extraction. 
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Figure 3. — System integration diagram. 

Front-end 

The front-end is the raw data collector and signal pre -processor. The B&K Pulse system was selected as the 
front-end for our system. The Pulse system is a real-time multi-analysis platform for vibration and sound analyses. 
The order analyzer of the Pulse system was activated in the integrations. In order to improve the signal-to-noise 
ratio, a synchronous averaging technique was also employed. The 1st order synchronous vibration was the main 
focus in the system. Multiple orders were monitored and processed in real time. 

Data Management 

Data management was carried out by the PTA (Product Test Adviser). PTA is a special version of SigQC 
developed by Signalysis Inc. SigQC was originally designed for manufacturer quality assurance. It is a production 
line test that responds to the demand among manufacturers to automate a process for accepting or rejecting units on 
the assembly line based on measured noise and vibration data. SigQC provides basic statistics, criteria, strategies 
and pass/fail methods. The major functions of the PTA include Database Management (Database tree) and 
Production Line Interface (traffic control). 


MATLAB Interface 

The PTA provides basic statistics, criteria, strategies and pass/fail methods, but it is very difficult to use these 
basic functions to construct a diagnostic algorithm for crack detection. MATLAB codes for mass and crack 
imbalance detection functions were developed according to equation (9) or equation (10) for resonant synchronous 
vibrations. The codes were integrated into the crack detection system by utilizing the interface provided in the 
SigQC. 


IV. Small Rotor Test 


To verify the crack detection system integration and algorithms, a Bently Nevada RK 4 rotor kit was used as a 
test bench. 


A. Setup 

The test rig diagram and instalment setup are shown in figures 4 and 5, respectively. A single rotor disk with 
two bearings was used in the setup. The rotor disk has a weight of 810 g and balance slots at a radius of 30.5 mm in 
22.5° increments around the circumference (see fig. 6). The rotor shaft has a diameter of 10 mm and the bearing 
span is 400 mm with an additional 200 mm overhang. With such a setup, the fundamental frequency of the system is 
at 2000 rpm. Modal tests indicate that the second bending resonance is at about 10,000 rpm. 
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Figure 4. — Small rotor rig test diagram. 



Figure 5. — Small rotor test kit. 



Figure 6. — Rotor kit disk with balance slots. 
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Proximity probes were used as the signal pick-ups. Sensor #1 (vibl) was in the vertical direction, while 
sensor #2 (vib2) was in the horizontal direction. A typical synchronous response of the system is shown in figures 7 
and 8, respectively. 

Nylon set screws, instead of steel ones, were used as imbalance weights to improve the imbalance resolution. 

As a result, the smallest imbalance unit was 3.416 g-mm. The disk 0° position was aligned with the center of the 
notch on the rotor shaft. A mechanical 1/rev signal (shaft notch detected by proximity probe) was used as the trigger. 
The rising edge of the signal was used in the trigging, which turned out to be the trailing edge of the shaft notch. 
Later a laser calibration indicated that the actual triggering point corresponds to approximately the 30° rotor disk 
circumferential position. 




Figure 7. — Typical synchronous response, vertical. 




Figure 8. — Typical synchronous response, horizontal. 
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B. Tests 


The initial system was balanced to the calibration system capacity. To verify the system integration and 
algorithms, imbalances were added at two circumferential positions, nominally 0° and 90° on the rotor disk, one at a 
time. The 1/rev signal comes from a notch on the shaft near the 0° mark. The actual notch is a filing-off of the shaft, 
which covers approximately 60° of circumference. According to the calibration, the imbalance correspond to 
approximately 30° and 120° with respect to the 1/rev trigger. At each circumferential location, 4 different imbalance 
weights were added: 1, 2, 3, and 4 set screws, which correspond to 3.416, 6.832, 10.248, and 13.664 g-mm 
imbalances. 

The motor speed was controlled such that the synchronous data acquisition between 350 to 8050 rpm with 
50 rpm step was guaranteed. Ten runs were recorded at the balanced state and a baseline was constructed based on 
these balanced inns. For each imbalance case, five repeated runs were recorded for further analysis. 

C. Data Reduction 

In the imbalanced runs, the baseline data were removed from each real time data acquisition by complex 
subtractions. The complex residue data were then fed in to the MATLAB based parameter identification algorithm 
through an iterative scheme. Both the amplitude and the phase of the imbalance were identified. Since the rotor shaft 
first bending was the dominant mode, a single-degree-freedom method was used, though multi-degree-freedom 
methods are also available in the codes. 

Typical parameter identification results are shown in figure 9. As a result of the high signal-to-noise-ratio near 
resonance, good curve-fitting results are usually achieved in both the amplitude and the phase. The fitting parameter 
is a complex number (fig. 10); its amplitude is proportional to the additional imbalance relative to the baseline, 
while its phase indicates the circumferential location of the additional imbalance with respect to the trigger. 




Shaft Speed, rpm 

Figure 9. — SDOF curve-fitting on the complex residue. 
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Figure 11. — Resonance-based identification-amplitude, vertical sensor. 


The imbalance amplitude and phase results from the vertical sensor curve fitting are shown in figures 1 1 and 12. 
It is seen that the amplitude identifications are approximately linearly proportional to the added mass imbalance, 
while the phase identifications are very consistent at the same circumferential location. The standard deviation of 
phase identifications is 4.45° at “0° mass imbalance” and 2.95° at “90° mass imbalance”. For the first 10 baseline 
runs, the amplitudes identified are small and the phases tend to be random, as expected. 

With known modal mass, the exact modal imbalance can be estimated without difficulty. For quick comparison, 
the identified parameters are scaled with respect to the mass imbalance at “0°” and “90°” at the smallest imbalance 
amount: 3.416 g-mm. The averaged scaled amplitudes are shown in figure 13. The percent change of the scaled 
amplitudes from that of the exact values is shown figure 14. For both circumferential locations, the percent changes 
are within 1 1 percent. 
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Similar data processing was carried out on the horizontal sensor. The results are shown in figure 15 to figure 18. 
Similar results are achieved. The percent changes are less than 9 percent in this sensor. 
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CD 

-Q 


200 


O) 

CD 

"O 

150 

c : 
o 

03 

100 

O 

O 

_l 

50 

c 

CD 

j— 

0 

£ 

E 

=3 

O 

-50 

i— 

b 

-100 


o e opo-euooeQGeeQOQQo 


GoeeoQoO^eooeeQeoec ) 


oo 


-oo- 


10 


20 


30 


40 


50 


Cycles 

Figure 16. — Resonance-based identification-phase, horizontal sensor. 


NASA/TM— 2006-214335 


11 


CD 

O 

c 

_CD 

CD 

-Q 

E 


TD 

Q) 

CD 

O 

CO 


E 

E 

■ 

CD 

<D 

"O 

13 


CL 

E 

CD 



Imbalance 


Figure 17. — Scaled amplitude factor, horizontal sensor. 


20 % 

15% 

5 10 % 

£ 5% 

a) 

f 0% 

£ -5% 

o 

® -10% 

-1 5% 

- 20 % 

3.4g-mm 6.8g-mm 13.7g-mm 20.5g-mm 

Imbalance 

Figure 18. — Comparison of amplitude to baseline, horizontal sensor. 


□ 0 deg 

□ 90 deg 



V. Engine Test 

A seeded fault engine test was also performed to demonstrate various crack detection strategies on an Allison 
TF41-A1B engine (ref. 8) as shown in figure 19. The GE system used vibration sensors (accelerometers) mounted 
on the engine case. The detection algorithm introduced in section 2 was implemented and integrated into the system. 
The advantage of this system is that the sensors used either already exist in an engine system or can be easily 
installed. 

The size and location for the initial fault seed was recommended by a stress analysis. It was determined that the 
location chosen was at a high stress area and the size and location of the slot was conducive to crack propagation 
(refer to fig. 20). 

The initial crack was a flaw of 0.005 in wide by 0.100 in radius half penny notch into the surface by an 
electrodischarge machining (EDM) process. The seeded fault was pre-cracked at a spin pit. The initial fault on the 
fan disk was propagated from each tip of the seeded crack about 100 and 105 mils, respectively, before installing in 
the engine for test (fig. 21) (ref. 9). 
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Figure 19. — Allison TF41-A1B engine on test stand. 



Figure 20. — Pre-cracked fan disk. 



Figure 21. — Pre-cracked seeded fault. 
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A. Test Setup 


The GE system for this engine test used four accelerometers as sensors. The list of sensors is as following: 
sensor #1: PCB 353B32, SN 83106, sensitivity 51.0mv/g; sensor #2: PCB353B32, SN 83107, sensitivity 51.1mv/g; 
sensor #3: PCB 353B04, SN 81329, sensitivity 10.48mv/g; sensor #4: PCB 353B04, SN 81330, sensitivity 
10.33mv/g. 

Sensors #1 and 2 were in the same plane as viewed from a cross section of the engine, and 90° apart. Likewise, 
sensors #3 and 4 were in the same plane and 90° apart. The sensors were mounted on the 1st stage stator case 
(fig. 22), where the vibrations from engine rotating parts can be best transmitted to the stationary engine case. 
Sensors #1 and 2 were used as main sensors. Sensors #3 and 4 were backups, in case of any overloading or other 
potential damages during the test. Sensor signals were conditioned at the test pad. The signal conditioners were 
protected by a steel cabin to prevent potential damage by an engine burst. The conditioned sensor signals were then 
transmitted to the control room for further on-line processing (fig. 23). A 1/rev tachometer signal and the 
acceleration signals were also recorded by a TEAC 8-channel digital tape recorder for future off-line analysis as 
needed (fig. 23). 



Figure 22. — Sensor locations. 



Figure 23. — Engine test system setup. 
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B. Data Processing 


The cycle profile used for the complete engine test was defined as: 27 short cycles plus one long cycle within 
half an hour. The short and long cycles are defined as: 

• The longer cycle consisted of a 30 sec dwell at idle, 75 sec ramp up to the maximum speed (8000 rpm), 

30 sec dwell at the maximum speed, and a 15 sec coast down to idle (approximately 2500 rpm) for a total 
cycle time of 150 sec. 

• The shorter cycle consisted of a 15 sec dwell at the idle, a 15 sec ramp up to the maximum speed 
(8000 rpm), a 15 sec dwell at the maximum speed, and a 15 sec coast down to the idle (2500 rpm) for a 
total cycle time of 60 sec. 

After the first day test, all the detection systems participating in the testing found no significant crack 
propagation. To increase the crack growth rate, the team decided to increase the maximum speed to 8900 rpm. 

Two weeks of testing finished with more than 4470 engine powered cycles without engine burst. The eddy- 
current inspection of the crack revealed a crack length of 0.130 ±0.010 in beyond the EDM notch on the outer 
diameter side of the EDM notch and 0.120 ±0.010 in beyond the inner diameter side of the EDM notch on the disk 
(fig. 24). This implies an approximately 25 mils radius averaged growth after two weeks’ testing. 

Figure 25 shows an averaged synchronous vibration response amplitude for sensor #1. El ere, the most 
significant resonant peak is located at about 5000 rpm. Thus, the MDOF system response model as described from 
equations (7) through (10) was used. Just the portion near the most significant resonant peak (5000 rpm) was used 
and a curve-fitting was carried out. As shown in figure 26, where the top portion of the figure is displaying the 
amplitude results while the bottom portion is the phase results, the red circles are the raw data from measurement 
and the blue line the curve-fitting results. A good fitting can be achieved by using equation (10) and the proposed 
iteration method. 



Figure 24. — Disk crack propagation after two weeks of testing. 
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Figure 25. — Averaged synchronous vibration response - sensor #1. 
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Figure 26. — Curve fitting results. 

The identified mass imbalance term and the crack-induced term are shown in figure 27. There are some cycles 
where the identified results are abnormal, largely due to imperfect data acquisition and diverging identification. 
Imperfect data acquisition usually happened after the starting up and shutting down of the engine. 

To better reveal the trend of the imbalance identifications, the abnormal points are removed to form the 
“cleaned data”. Then, a wavelet transformation technique was applied to the “clean data” to extract the long-term 
trend as well as day-to-day variational trends. A level 12 Daubechies wavelet was used in the transform. The overall 
trend and the day-to-day trend, together with the “cleaned data” are shown in figure 28. In figure 29, the overall 
trend and day-to-day variational trends are extracted and tagged with date marks. 

As revealed in detail in figure 29, only minor increases happened in the over-all trend during the 2 weeks’ 
testing. The day-to-day trend indicated the variation within a day, with a pattern of low-to-high-to-low. This is very 
likely due to dramatic daily temperature changes in the desert location where the tests were conducted. 

Similar analyses on sensor #2 can be found in figures 30 through 32, where the overall trend is slightly 
decreased and the day-to-day trend exhibits the same variations in the trends from sensor #2. 
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Figure 27. — Sensor #1 identification results. 
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Figure 28. — Sensor #1 abnormal removed data and trends. 
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Figure 29. — Sensor #1 overall trend and day-to-day trend. 
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Figure 30. — Sensor #2 identification results. 
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Figure 31. — Sensor #2 abnormal removed data and trends. 
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Figure 32. — Sensor #2 overall trend and day-to-day trend. 
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VI. Conclusions 


A method for rotor imbalance identification using rotor resonant synchronous vibration was presented. The 
imbalance is identified as a complex number through an iterative least-mean-square-based curve fitting technique. 
The magnitude of the identified complex number is the magnitude of the modal imbalance, while the phase of it is 
associated with imbalance location in the circumferential direction. 

A complete data acquisition, data reduction, and data management system was integrated using state-of-the-art 
hardware and software. A MATLAB script was developed and interfaced with the system for real-time imbalance 
detections. A small rotor test rig was successfully used to validate the system integration and to examine the 
proposed algorithms. 

The developed system was applied to a real engine test with embedded disk crack. The engine disk crack did 
not propagate into engine burst in the planned two weeks’ testing. Though the crack did propagate for 45 mils, the 
crack growth was not large enough to be detected by the current method. The system data reduction using a wavelet 
transform did reveal minor day-to-day variations, which are very likely due to large temperature variation in the 
testing site. 

Because of high signal-to-noise-ratio at resonance, the resonance-based algorithm is expected to have higher 
sensitivity than other vibration based method. It can be devised to use in health monitoring system for rotating 
machinery. It can also serve as an alternative for rotating machinery balancing. Further investigation of this 
technology is recommended. 
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